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ABSTRACT 

We present the results of a study of a sample of 375 Extremely Red Galaxies (ERGs) in the Phoenix 
Deep Survey, 273 of which constitute a subsample which is 80% complete to K s — 18.5 over an area 
of 1160 arcmin 2 . The angular correlation function for ERGs is estimated, and the association of ERGs 
with faint radio sources explored. We find tentative evidence that ERGs and faint radio sources are 
associated at z > 0.5. A new overdensity-mapping algorithm has been used to characterize the ERG 
distribution, and identify a number of cluster candidates, including a likely cluster containing ERGs at 
0.5 < z < 1. Our algorithm is also used in an attempt to probe the environments in which faint radio 
sources and ERGs are associated. We find limited evidence that the I — K s > 4 criterion is more efficient 
than R — K s > 5 at selecting dusty star-forming galaxies, rather than passively evolving ERGs. 

Subject headings: galaxies: clusters: general - galaxies: high-redshift - galaxies: structure - infrared: 
galaxies - surveys 



1. INTRODUCTION 

Extremely red objects (EROs) comprise stars, stellar 
remnants and galaxies. EROs that are not galaxies in- 
clude bro wn dwarfs and p rotostars, and have been stud- 
ied (e.g. iFan et all [20001 ) at infrared K magnitudes < 
15. At fainter magnitudes, most EROs are galaxies. Use 
of the term "ERO" to describe both classes of objects is 
widespread, but to avoid confusion, we refer hereafter to 
the ERO galaxies as extremely red galaxies (ERGs) . They 
have very red colours (R — K s > 5, I — K s > 4), and are 
thou ght to in clude members of two distinct popu lations 
(e.g. iCimattil [2003L iMcCarthvl [200l iDaddi et alJl200l : 
dusty star-forming galaxies (including some active galac- 
tic nuclei, AGNs) an d evolved ellipticals at z > 1. Fol- 
lowing the example of iRoche et al.l (|2002f) . we will refer to 
the dusty star-forming ERGs as "dsfERGs", and to the 
passively evolving ERGs as "pERGs". The relative mix 
of the two populations, which is a function of magnitude, 
colour and redshift, has not been well constraine d in the 
literature, although progress has been made (e.g. ICimattil 
120031 iMoustakas et alj|200l . The red colours in dsfERGs 
are caused by dust absorbing shorter wavelength light and 
in pERGs by old stellar populations. 

Different models of galaxy evolution — monolithic col- 
lapse with passive luminosity evolution (PLE) versus hi- 



erarchical merging in a cold dark matter (CDM) cosmol- 
ogy — predict different properties for galaxy populations. 
These involve differences in the formation scenario for el- 
lipticals, and in the evolution of the large-scale structure 
of the universe. The hierarchical merging model predicts 
that the number of ellipticals should increase with time, 
as they form through the merging of spiral or irregular 
s ystems. 

lEisenhardt et alj (|2000f ) found that a high fraction of 
red galaxies, corresponding to redshifts z > 1, was at 
least as consistent with P LE as with CDM merger models. 
Scodeggio & Silva (|2000[) found that the space density of 
ERGs was consistent with no change in the volume density 
of ellipticals from z = to 1.5. However, Rodighiero et 
al. (2001) found evidence for CDM merger models in the 
form of a decrease in the density of E/SO galaxies between 
z — 1 and 1.5, with a corresponding drop in the density 
of ERGs. A characteristic limitation of these early ERG 
surv eys was large field-t o-field variation s in ERG density 
(e.g. iBarger et al~l 119991 compared with iMcCracken et all 
2000). This is largely a result of cosmic variance, but is 
also rela ted to survey sen sitivity and how the ERG class is 
defined ([McCarthy! l2004h . ERG density is a strong func- 
tion of both the survey depth and the colour threshold 
used to define the ERG population. 
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ERGs include galaxies that for med at high redshi ft and 
evolved passively since then (e.g. iFirth et al.ll2002t ). For 
example, Spinrad et al. (|1997l ) studied a pERG at z = 1.55, 
estimating that it had a formation redshift Zf > 5 and 
identifying it as possibly the oldest galaxy known at z > 1. 
The detection of highl y evolved ellipticals at 1 < z < 2 by 
Bcm'tcz et al. ((l999), with a corresponding lack of lu- 
minous blue ellipticals, also implied a formation redshift 
for ellipticals of Zf > 5. From ERGs in the K20 survey, 
Cimatti et al. (|2003[ ) deduce a minimum formation red- 
shift of about 2 for pERGs at z ~ 1, and a significant 
range for the formation redshift of about 2.2 to 4. 

Even though the evolution of the star-formation rate 
(SFR) in different environ ments has not been comprehen- 
sively constrained fe.g. iPoggianti et ail I2006D . it is clear 
that the lowest SFRs in the nearby universe, found in high 
mass elliptical galaxies, were the s ites of the highest SFRs 
in the distant past. Smail et al. (|l999l ) suggest that the 
most extreme ERGs are mainly of the dusty star-forming 
variety, and that studying the dsfERG population may tell 
us much about the nature of star formation in the early 
universe. 

Measurement of the densities of pERGs may therefore 
be used to test models of galaxy evolution in terms of the 
mass assembly of galaxies, with the clustering properties 
of both pERGs and dsfERGs being used to test the clus- 
tering predictions of these scenarios. The importance of 
ERGs for such tests is emphasized by Georgakakis et al. 
(2006), who find that ERGs contribute almost half the 
stellar mass density of the universe at z rj 1. Further- 
more, both ERG classes may provide valuable information 
about the formation of the earliest galaxies and clusters. 
Faint radio so urces are also consid ered to identify clus- 
ters (e.g. Best 120001 Buttery et al. I2003D . but it has not 
been conclusively shown that ERGs and faint radio sources 
trace similar structures at high redshift (Georgakakis et al. 
120051) . 

In this paper, we present the results of a moderately 
deep near- infrared K s survey (Section 2) covering « 1160 
arcmin 2 in the Phoenix Deep Survey region. Our survey 
covers a larger area than the -ftT s -band survey detailed in 
Sullivan et al. (|2004h but to a shallower depth. Our corre- 
sponding ERG sample, defined in Section 3, is about the 
same numer ical size as that studied by Georgakakis et al. 
( 2005], I2006D . Our ERG sample is complemented by deep 
radio data, which we use to examine the association of 
ERGs and radio sources. Deep multiwavelength coverage 
allows us to investigate differences in clustering properties 
which may be introduced by different ERG selection cri- 
teria. We examine the environments and clustering of our 
sample in Sections 4 and 5, focussing on a possible cluster 
at 0.5 < z < 1 in Section 6. Where necessary in this work, 
we have adopted a concordance cosmological model with 
H = 70 km s- 1 Mpc" 1 , Q M = 0.3 and 0\ = 0.7. All 
magnitudes are in the Vega system. 

2. THE PHOENIX DEEP SURVEY 



The Phoenix Deep SurvejEl (PDS) is a multiwavelength 
survey aimed at studying the nature and evolution of faint 
radio sources. The radio data, which reaches a minimum 
1 a noise level of 12 /xJy at its most sensitive, defines 
the survey area of 4.56 deg 2 in the southern constel lation 
Phoenix. T he radio observations (Hopkins et al. .1998, 
119991 I2003D were carried out between 1994 and 2001 at 
the Australia Telescope Compact Array (ATCA). 

2.1. Survey Multiwavelength Data Set 

Deep multicolour images are available for part of the 
PDS in U, B, V, R and / bands. BVRI data were obtained 
in 2001 using the Wide Field Imager (WFI) at the Anglo- 
Australian Telescope (AAT), Siding Spring Observatory. 
U band data were collected at the Cerro Tololo Inter- 
American Observatory (CTIO) in 2003, while deep near- 
infrared data (K s < 20) are acquired for a 180 arcmin 2 
region using the Son of ISAAC (Sofl) instrument at the 
European Southern Observatory's (ESO) New Technol- 
ogy Telescope (NTT) at La Silla. Details of the optical 
and n ear-infrared observations are given in Sullivan et al. 
(2004J). The R— and I— band catalogs have 5 a limiting 
magn itudes of 24.61 and 24.67 respectively, using SExtrac- 
tor's (|Bertin fc Arnoutslll996l ) MAG_AUTO photometry. 

2.2. Extended Near-Infrared Dataset 

In addition to the Sofl data, near-infrared images (K s < 
18.5) covering an area of 1160 arcmin 2 of the PDS were 
obtained with the Infrared Side Port Imager (ISPI) cam- 
era at CTIO on 2004 September 25. The conditions were 
photometric, with seeing ranging from 0'.'8 to l'/2. We ob- 
served 11 ISPI pointings using a random dither pattern, 
nine with a total exposure time of 30 minutes, and the re- 
maining two with exposure times of 21 and 15 minutes. 
Image processing was performed using standard NOAO 
IRAB3 routines. Dark frames were subtracted from each 
image, as were running sky flats, constructed from medi- 
ans of up to 10 object frames. The IRAF tasks mscgetcat 
and msccmatch were used to accurately register the images 
for co-addition. The final astrometry w as done by refer- 
ence to 2MASS3 dSkrutskie et al.l[2CJ06T) source positions. 
Residual offsets with respect to 2MASS positions, and to 
fainter PDS sources in i?-band images, were fitted by a 
Gaussian with a — 0'.'3. 

We employed SExtractor on each of the eleven if s -band 
images using their respective exposure maps as weight im- 
ages, with a detection threshold of 1.5 times the noise level. 
The low signal-to-noise regions on the edges of the co- 
added images, a result of the dithering of individual expo- 
sures, are excluded with simple right ascension and decli- 
nation cuts. The overlap between the images was such that 
we were able to maintain contiguity throughout the area 
of our if s -band survey. The final catalogue contains 4834 
objects. Within the area of our catalogue, there are 230 
sources, mostly stars, which also have 2MASS identifica- 
tions. For each individual image, we tied our photometry 
to 2MASS by adjusting the magnitude zeropoints in our 



1 See also http://www.physics.usyd.edu.au/~ahopkins/phoenix/ 

2 IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Association of Universities for Research in 
Astronomy, Inc., under cooperative agreement with the National Science Foundation. 

3 This publication makes use of data products from the Two Micron All Sky Survey, 2MASS, which is a joint project of the University of 
Massachusetts and the Infrared Processing and Analysis Center /California Institute of Technology, funded by the National Aeronautics and 
Space Administration and the National Science Foundation. 
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SExtractor input to give the tightest correlation between 
ISPI and 2MASS magnitudes (see F igure CD). 

Following Schmidt et al. (|2006t) . we performed star- 
galaxy separation by comparing SExtractor's fixed aper- 
ture (3'.' 5 diameter) magnitude with its MAG.AUTO mag- 
nitude, after which we used MAG_AUTO for all subse- 
quent analysis. Our _ftT s -band galaxy counts shown in Fig- 
ure [2] are consistent with other recent measurements. The 
completeness of the catalogue was estimated to be 80% to 
K s = 18.5 fr om comparison with source counts from Kong 
et al. (|2006D . 

3. THE ERG SAMPLE 

3.1. ERG Selection 

We matched our ne ar-infrared catalogue with the exist- 
ing optical catalogue (jSullivan et al.1 [20041 ) , counting ob- 
jects within 2" as matches. We constructed two samples 
of ERGs based on the two-colour criteria: R — K s > 5 and 
I — K s > 4. This resulted in a sample of 375 R — K s - 
selected ERGs, 273 of which are brighter than K s = 18.5. 
Similarly, we derived a sample of 346 I— if s -selected ERGs, 
256 of which are brighter than K s =18.5. 301 ERGs are 
common to both samples, 228 of which are brighter than 
K s =18.5. Figure [3] shows a comparison of the differen- 
tial R — K s ERG counts with those from previous studies. 
ERGs with more extreme colours close to the K s -han& 
detection limit may not be detected in i?-band, and will 
therefore contribute to the incompleteness seen in Figure 

El 

3.2. ERG Correlation Function 

To characterize the clustering properties of this ERG 
sample, we use the two-point angular correlation func- 
tion. We fo llow the prescription detailed by G eorgakakis 
et al. (|2005[ ). using the Landy & Szalay (|1993f ) estimator 
for calculating w(0), with uncertainties that are assumed 
to be Poissonian. We also apply an integral const raint 
(C), again as detailed by Georgakakis et al. (|2005[ ), as- 
suming the correlation function slope 5 = 0.8, which gives 
C = 5.42 for our survey area. 

In Figure 2] we present the correlation functions for the 
ERGs from our data. We find significant clustering for 
limiting K s magnitudes between 18 and 18.5 (Table [1]). 
The clustering amplitudes for our ERG sample are con- 
sistent with other recent measurements (Figure [5J , and 
are much greater than the amplitudes which have been 
measured for all if s -band galaxies t o the same limiting 
magnitudes (e.g. iDaddi et al.ll2000T) . There appears to 
be a difference between the clustering properties of the 
ERGs selected by our two different criteria. In Table [Q 
at each limiting magnitude, the clustering amplitude and 
confidence level of the clustering signal for ERGs selected 
using the R — K s criterion is greater than that found using 
I — K s selection. This is also seen in F igure [H supporting 
the suggestion (Vaisanen & Johansson 2004) that I — K s 
selection may be more sensitive to the dsfERG population, 
as these galaxies are expected to cluster less strongly than 
the pERGs. We stress that differences in the redshift dis- 
tributions of ERGs selected by the different colour criteria 
are likely to affect their measured clustering prop erties . 
Observed-frame colour tracks given by McCarthy (|2004h 



show that the R — K s criterion should, in principle, be 
more efficient at selecting passively evolving ERGs. 

4. ERG ENVIRONMENTS 

Studying the association of ERGs with i4T s -band galaxies 
and faint radio sources is valuable in understanding their 
nature. For this analysis we have used the ERG sample 
selected with R — K s > 5. 

Our estimation of the cross-correlation between K s - 
band galaxies and ERGs reveals no significant association. 
Given the depth of our survey, in which most of the K s - 
band galaxies lie at z <C 1, this is not unexpected, since 
ER Gs are found predominantly at z « 1. Georgakakis et 
al. (|2005| ) find that ERGs are only associated with over- 
densities of i^-band galaxies at z > 1. 

To compare the ERG and faint radio source populations, 
we first make use of th e publicly available code hyperz 
(Bolzonc lla et al.ll2000t ) to find photometric redshifts for 
all galaxies detected in at least three of the UBVRIK S 
bands. In Figure[S]we have shown the photometric redshift 
distribution for R — K s ERGs only. We then match our 
radio catalogue with the galaxies (including non-ERGs) 
for which we have photometric redshifts, dividing the ra- 
dio sources into two broad redshift ranges. Faint radio 
sources individually associated with ERGs are excluded, 
and cross-co rrelati on functions are estimated (e.g. Geor- 
gakakis et al. 120051) to compare the two radio samples with 
the ERG catalogue. Figure in which bootst rap errors 
are shown (Barrow, Bhavsar & Sonoda Il984j ). demon- 
strates that while there is no significant correlation be- 
tween the positions of ERGs and faint radio sources with 
counterparts at z p hot < 0.4, there is such a correlation at 
the 2.3ct confidence level for the radio sources with coun- 
terparts at £phot > 0.5. Georgakakis et al. (|2005f ) also 
detect such a correlation at the w 2 a confidence level. 

Our work supports the conclusion that at higher red- 
shifts [z > 0.5), ERGs and radio sources are associated. 
This is perhaps not surprising since the redshi ft dist ribu- 
tions o f both samples peak at z « 1 (Condon |l989l Mc- 
Carthy |2003). Since wc have eliminated faint radio sources 
individually associated with ERGs, this cannot simply be 
the result of one-to-one associations. If ERGs and radio 
sources trace galaxy overdensities at higher redshift, we 
should be able to detect strong association between the 
two populations in high redshift cluster environments. In 
the next section we search for this effect. 

5. CLUSTER CANDIDATES 

Rich clusters at high redshifts contain massive, red el- 
liptical galaxies at their cores, which can show up as ERG 
overdensities a few ar cminutes acro ss in E RG surveys (e.g. 
St anford et al. Il997l Roche et al. l2002t Diaz-Sanchez et 
al. H007). Studies of such cluster samples may be used 
to investigate how the highest density regions of the uni- 
verse have evolved, and possibly place constraints on their 
formation. 

Here we identify cluster candidates in two different ways. 
First, we use a simple counts- in-cells criterion, and sec- 
ond, we use an overdensity-mapping algorithm. The for- 
mer technique gives us a set of likely candidates, and the 
latter technique allows us to treat cluster environments 
statistically. 
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5.1. Counts- In- Cells 

For the entire area of our survey we have found local 
counts of ERGs within circles of radius 50 arcseconds, 
and identified six locations (Figure [8|) which contain six or 
more ERGs wi thin a circle of radius 50 arcseconds (see also 
Blake & Wall |200"3 ). Given a uniform random field with 
the same mean density as our ERG catalogue, the prob- 
ability of a particular location hosting such an overden- 
sity is 3.8 x 10 -4 . This criterion selects locations hosting 
the highest densities and is intended to identify overden- 
sities with the highe st pro bability of physical significance. 
Diaz-Sanchez et al. f|2007T ) report one such overdensity as 
a cluster at z « 1. Figure [9] shows the distribution of faint 
radio sources in these cluster environments, for the same 
redshift ranges used in Figure [7] One might expect from 
the statistically significant result for z > 0.5 in Figure 
[7J that there should also be an excess of associated radio 
sources in the regions with ERG overdensities. However, 
we find only a 1.05 a excess of faint radio sources in these 
putative clusters. This may be because the ERGs in these 
high density regions are pERGs which have fewer radio 
counterparts since they have lower star-formation rates. 

5.2. Overdensity- Mapping Algorithm 

We now describe an algorithm we have developed to 
identify overdensities of ERGs within our survey area, be- 
ginning with a general description before explaining the 
procedure in detail. For any given location, we compare 
the local galaxy count with the distribution of counts ex- 
pected from a large set of uniform random catalogues, each 
with the same mean density as the ERG catalogue. This 
comparison allows us to calculate the probability that the 
local area contains a greater number of galaxies than that 
expected from a random distribution. Repeating this pro- 
cess for a range of specific radii, we form a probability func- 
tion which is then associated with the location. Similarly, 
a probability function is assigned to each of a set of uni- 
formly spaced locations within our survey. The optimal 
range and resolution of the radii and locations sampled 
may be found empirically. 

Each probability function is used to estimate the proba- 
bility of the count within each radius enclosing its location 
occurring as the result of a uniform random distribution. 
If the probability of obtaining the local count from a ran- 
dom catalogue is high, the value of the probability function 
will be close to a half, since the number of random cata- 
logues with lower counts will be similar to the number with 
higher counts. In the presence of a genuine overdensity, 
the probability function will rise above 0.5 and continue 
to ascend with increasing radius (since the probability of 
a large overdensity is lower than the probability of a small 
one) until the "edge" of the overdensity has been reached. 
Beyond this edge, the enclosed density will fall back to 
the background density, and the probability function will 
approach 0.5. Hence, we expect each probability function 
to be sensitive to both the significance and scale of any 
nearby overdensities, characterised by its maximum and 
the radius at which it occurs. 

For overdensities centred on particular catalogue ob- 
jects, this maximum occurs at a radius of zero (which thus 
encloses an infinite density) , but we negate this by specify- 
ing a minimum overdensity membership of two. We define 



the radius at which the maximum occurs as the overden- 
sity radius, and the overdensity probability as the value 
which the probability function takes at the overdensity ra- 
dius. Our technique uses probability to compare between 
scales in order to select the most likely scale at which each 
particular location hosts an overdensity. 
To achieve this, we follow these steps. 

• The survey area is covered with a uniform grid of 
test points at intervals of 20 arcseconds along both 
right ascension and declination. 

• We populate the survey area with a large number 
of random sets, each containing the same average 
density of points as the ERG catalogue that we are 
examining. 

• We count the number of ERGs and random objects 
(in each random set) as a function of distance from 
each test point, in steps of 10 arcseconds from 20 
to 120 arcseconds. 

• The probability of overdensity is then the fraction 
of random sets for which the number of ERG counts 
exceeds the number of random counts. ERG counts 
consistent with uniform random distribution will 
therefore produce overdensity probabilities of 0.5. 

In overdense regions, the overdensity probability directly 
quantifies the probability that the local environment is as- 
sociated with an overdensity, and the overdensity radius 
directly characterizes the scale size of this structure. Each 
of the grid points which we test within our catalogue will 
have an associated overdensity probability and overdensity 
radius. By confining our random sets to the area of the 
catalogue under inspection, we account for edge effects. 

The probability of overdensity is intimately connected 
with the overdensity radius. We are effectively sensitive to 
multi-scale structure throughout the survey area, meaning 
that we will detect both high density small-scale features 
(which might be overlooked in the case of simple large- 
scale smoothing), and extended regions of intermediate 
density (which might be overlooked in the case of simple 
small-scale smoothing). Our algorithm is therefore sensi- 
tive to overdensities spanning a wide range of scales. In 
particular it has the potential to simultaneously identify 
filamentary structures as well as clusters, an aspect that 
we will explore in detail in future work. 

Figure [10] shows the result of the overdensity-mapping 
algorithm applied to our data. It is clear that the regions 
identified by our algorithm as having a high overdensity 
probability do indeed correspond to visual overdensities of 
ERGs, and similarly that regions with no ERGs, or iso- 
lated ERGs, have low probability. 

5.3. Faint Radio Sources and the Overdensity Map 

In Figure [TT] we visually compare the overdensity proba- 
bility map in Figure [Tm with the distribution of faint radio 
sources in the two z p hot ranges previously defined. The 
redshift dependence shown in Figure [Jj is not apparent in 
Figure [TTJ We stress, however, that since our algorithm 
has smoothed the predominantly z ~ 1 ERG population, 
the overdensity map will be biased towards structures at 
this redshift. 
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We have matched each of the faint radio sources with 
the pixel (from our overdensity map) closest to its position, 
and counted the frequency with which faint radio sources 
are matched with certain pixel values. Effectively, what 
we are doing here is counting the number of faint radio 
sources which lie along different probability contours of 
our overdensity map. The result of this analysis is shown 
in Figure [TSJ which suggests that faint radio sources with 
optical counterparts with z p hot > 0.5 are indeed found 
in high density regions, with a 2er excess for overdensity 
probabilities between 0.9 and 1. This is not seen in the 
•Zphot < 0.4 sample. While the statistics are poor, and our 
overdensity map is strongly biased towards higher redshifts 
(as shown in Figure [6]), Figure [T2l tentatively suggests that 
clusters at high redshift are more likely to host faint radio 
sources in their central regions. 

The lack of a clear association of ERGs and radio sources 
with optical counterparts with z p hot > 0.5 evident in Fig- 
ure [TT] and the result from Figure [9] suggests that the cor- 
relation between ERGs and moderately high redshift radio 
sources implied by Figure [7] does not occur in regions with 
the highest ERG densities. Figure [7] shows association on 
scales which may not be noticeable when considering the 
scale of the entire field shown in Figure [Til an d cannot be 
used to draw conclusions about the characteristic density 
of the environment in which ERGs and faint radio sources 
are associated. Figure [TJ] illustrates, however, that our 
overdensity map provides a way of directly probing this 
characteristic density. 

6. A POSSIBLE CLUSTER AT 0.5 < Z < 1 

Of the cluster candidates identified via the counts-in- 
cells criterion (Figure [5]), we have singled out one in par- 
ticular as an interesting object for future study. 

-K^-band and i?-band images of this putative cluster are 
shown in Figure [TBI with ATCA 1.4 GHz contours over- 
laid. Astrometry for both images was refe renced to 8-10 
stars in the SuperCOSMOS (Hambly et al. 120011 ) i?-band 
image which has formal fit errors of 0.11 arcsec rms in 
each coordinate. The bright (K s = 16.3) central radio 
source PDF J011053.8-454339 has a flux density of S1.4 = 
106 /iJy. At this flux density, the relative proportions of 
starburst galaxies and AGNs in the faint radio popu lation 
are stro ngly m odel dependent (e.g. Seymour et al. 2004, 
Hopkins 2004), although recent observations (Seymour et 
al. 120081 ) indicate that star-forming galaxies at this flux 
density are more numerous than AGNs by a factor of 2 
to 3. Figure [TBI shows that the central radio source has a 
nearby companion radio source, PDF J011053. 4— 454351, 
associated with a faint i^-band object that falls below our 
catalogue threshold. 

Figure [14] shows there are six objects in this possible 
cluster with R — K s > 5.4. Since such extreme colours are 
caused primarily by redshift, the cluster is likely to lie at 
0.5 < z < 1 but not much greater, given the depth of our 
K s survey. If it lies at z = 1, the radio luminosity of the 
central ERG is L1.4 = 5 x 10 23 WHz -1 , making it a low- 
luminosity AGN. As the brightest galaxy in this field, it 
is also consistent with being the central dominant cluster 
galaxy. 

An 80 /zJy radio source to the west (PDF J011052. 0-454338) 
is also shown in Figure [TBI with contours suggestive of a 



head-tail morphology. Tailed radio gal axies a re often seen 
in galaxy cluste rs (e.g. Gomez et a l. Il997l Sakelliou & 
Merrifield 120001 Klamer et al. 120041) . but its R = 20 .4 
magnitude and clear separation from the red locus in 
Figure [JJ] point to it belonging to the foreground cluster 
evident in the i?-band image. 

Figure [T5] shows that the ERGs associated with this 
overdensity all fall within 22 arcseconds of the central 
ERG, corresponding to a projected radius of ~0.18 Mpc at 
z ~ 1 . This is less than 1/5 of an Abell radius (e.g. lAbelH 
1965). However, the radial extent of the cluster is likely 
to be underestimated because of our relatively bright K s - 
band limiting magnitude. The presence of radio sources 
within our empirically determined radial extent suggests 
that this estimate is consistent with the core radius of 167 
kpc found for nearby clusters, base d on their radio source 
distribution (Ledlow & Owen ll995l ). 

7. SUMMARY AND FUTURE WORK 

Within the Phoenix Deep Survey we have identified a 
sample of ERGs over a large area (1160 arcmin 2 ) that is 
80% complete to a limiting magnitude of K s — 18.5. The 
number counts and clustering properties of our ERG sam- 
ple are consistent with previous observations. Based on 
photometric redshifts we find evidence for the association 
of ERGs and faint radio sources with z p hot > 0.5, but 
not with Zphot < 0.4, consistent with earlier results on the 
evolution of star formation and AGN activity. 

We find weak evidence to suggest that R — iC s -selected 
ERGs are more strongly clustered than / — /sT s -selected 
ERGs. If confirmed with more extensive datasets, this 
would imply that R — K s -selected ERGs are more strongly 
associated with overdensitie s, sup porting the suggestion 
by Vaisanen & Johansson (|2004f ) that R — if s -selected 
ERGs contain a higher proportion of passively evolving 
sources, assuming that the two ERG populations have sim- 
ilar redshift distributions. 

The identification of the cluster candidates in this study 
shows that ERGs can be used to identify overdensities of 
potential physical significance, and our 2.3 a result for the 
association of ERGs and faint radio sources at z p hot > 0.5 
is evidence that both classes of objects trace overdense 
regions. As passive galaxies dominate in high density en- 
vironments, the absence of radio source overdensities in 
most of our clu ster ca ndidates supports the result s of Geor- 
gakakis et al. (|2006l ) and Simpson et al. ((2006) that ds- 
fERGs are more closely associated with the radio popu- 
lation than pERGs. As a result, faint radio sources are 
most efficient at identifying starburst galaxies in the infall 
regions of clusters. The probable 0.5 < z < 1 cluster iden- 
tified here will be explored in more detail with spectra and 
deeper imaging. 

We have introduced a method for quantifying overden- 
sities in a galaxy distribution that is sensitive to overden- 
sities on a broad range of scales, and demonstrated its 
efficiency in identifying cluster candidates. A faster im- 
plementation of this algorithm is being developed. 

ERG sample size is a strong function of limiting K s - 
band magnitude, and it is necessary to have complemen- 
tary deep optical R- and /-band data to detect ERGs with 
the most extreme colours. Large samples of ERGs are be- 
coming available with the next generation of large area 
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near- infrared surveys, such as UKIDSS (Lawrence et al. 
I2007T ) . Where these overlap with existing deep radio sur- 
veys, it will be possible to explore the association of the 
faint radio population with ERGs in unprecedented detail. 

We thank Helen Johnston, Scott Croom and Chris 
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Table 1 

Correlation Function Amplitudes A 







K s limiting magnitude 


Selection 


Properties 3, 


18.0 


18.25 18.5 


R-K s 


A (xlCT 3 ) 


28.7 ±4.1 


19.6 ±2.8 15.0 ±2.1 




N 


143 


202 273 


I-K s 


A (xl(T 3 ) 


26.9 ±4.4 


16.2 ±3.0 9.6 ±2.2 




N 


132 


188 256 



Uncertainties in A are calculated from 1 a errors in the fit- 
ted amplitudes. The number TV of objects brighter than each 
limiting magnitude is also shown. 
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Fig. 1. — A comparison of the 2MASS and ISPI magnitudes after photometric calibration, for objects with K s < 15. The points 
significantly below the line arise from extended or saturated objects, resulting in erroneous integrated fluxes. The error boxes indicate the 
typical uncertainties as a function of magnitude. 
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Fig . 2. — Differential galax y source counts, compared with counts from Georgakakis et al 
(1998) and Kong et al. 1 120061) . Differential star counts are also shown. The large error bar for the Georgakakis et al. point at K 



), Schmidt et al. (2006), Szokoly et al. 

12.25 



arises from a single object in that bin. The high counts at the bright end of our sample may arise from a small systematic overestimation of 
brightness by SExtractor at these levels. 
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Fig . 3. — Differential ERG number coun ts for the R — K a > 5- selected sample compared with Georgakakis et al. (2005), Daddi et al. 
(2000), FLAMEX/SDSS (Elston et al. 2006) and Kong et al. (2006). Our counts at the bright end are subject to small number statistics. 
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Fig. 4. — Correlation functions for the ERG samples as a function of limiting magnitude. Left: R — K a > 5 and Right: I — K s > 4. The 
value of <5 is fixed at 0.8 and the upper and lower lines signify ±1 cr uncertainties. The increased scatter in the I — i<" s -selected samples, 
introduced by a small number of objects that are not common to both samples, illustrates the difference between the clustering properties of 
the ERG samples selected by the two criteria. 
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Fig. 5. — Clustering amplitu des for the R — K a > 5 and / — K s > 4 samples compared with those found by Georgakakis et al. (2005), Daddi 
et al. <|2000f ) and Roche et al. l|2002T ) . For clarity, points at the same limiting magnitude have been offset horizontally by ±0.04 mag. At a 
limiting magnitude of K a = 18.5 our clustering amplitudes are significantly lower than Daddi et al. because of our survey incompleteness. 
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Fig. 6. — Photometric redshift distribution for the R — K 3 > 5 ERG sample. 
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Fig. 7. — Cross-correlation functions with bootstrap errors for ERGs around faint radio sources with optical counterparts assigned photo- 
metric redshifts (not shown in Figure[6]l. Left: 97 radio sources associated with galaxies at z p hot < 0.4 and Right: 98 sources with z p ^ ot > 0.5. 
As discussed, our exclusion of faint radio sources individually associated with ERGs ensures that the small-scale feature is not the result of 
one-to-one associations. 
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Fig. 8. — Our ERG field showing the circled locations of the cluster candidates based on our counts-in-cells criterion (see section 5.1). The 
cross indicates the probable 0.5 < z < 1 cluster discussed in section 6. 
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Fig. 9. — Density profiles for faint radio sources around the six positions (see Figure [8} which lie at the centres of overdensities selected via 
the counts-in-cells criterion. Left: z p ^ ot < 0.4 radio sources, and Right: z p i lot > 0.5 sources. For this measurement we have not eliminated 
objects that are both ERGs and faint radio sources. The horizontal lines in each panel correspond to the average radio source density. 
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Fig. 10. — Grayscale representation of our overdensity probability map for the Ft — K" s -selected ERG sample (see section 5.2). The black 
dots are the ERGs in this sample and the contours correspond to probabilities of 0.5, 0.625, 0.75, 0.875 and 1. There is clear evidence of 
filamentary structure in the high density islands. 
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Fig. 12. — Surface density of radio sources as a function of overdensity probability. Left: Sources with z p hot < 0.4. Right: Sources with 
z phot > 0-5. The error bars are Poissonian and the horizontal lines correspond to the average source density. 
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Fig. 13. — Images of a possible 0.5 < z < 1 cluster (see Section 6). Left: A" s -band; Right: R-band. Overlaid contours at 1.4 GHz from the 
Phoenix Deep Survey are at 36, 47, 61 and 79 /ijy/beam. The ATCA synthesised beam (12" X 6" in position angle 0°) is shown in the lower 
left hand corner of the Ks-band image. 
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Fig. 14. — Colour-magnitude diagram for the galaxies in the probable 0.5 < z < 1 cluster (Figures l8l and \13\ within 22" of the brightest 
alaxy. The brightest galaxy is marked with an asterisk, and the western tailed radio source host is marked with a circle. 
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Fig. 15. — The vicinity of the probable 0.5 < z < 1 cluster, showing the locations of cataloged i^ s -band galaxies (O)i faint radio sources 
(+) and ERGs (x). Where available, we have labeled each object with K a magnitude, R — K a colour and radio flux density. A circle of radius 
22 arcseconds is also shown. 



